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Abstract 

We study the magnetic properties oi S = 1/2 antiferromagnetic Heisenberg 



> 

o 

Q\ ' chains with inhomogeneity of interaction. Using a quantum Monte Carlo 

^^ ■ method and an exact diagonaUzation method, we study bond-impurity effect 

jrt ■ in the uniform S = 1/2 chain and also in the bond-alternating chain. Here 

'bond impurity' means a bond with strength different from those in the bulk 

. , or a defect in the alternating order. Local magnetic structures induced by 

^ ' bond impurities are investigated both in the ground state and at finite tem- 

O ■ peratures, calculating the local magnetization, the local susceptibility and the 

local field susceptibility. We also investigate the force acting between bond 

impurities and find the force generally attractive. 
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I. INTRODUCTION 

Quantum effects are very sensitive to the inhomogeneity of interaction. The magnetic 
structure induced by such inhomogeneity has much attracted us for its pecuhar properties. 
The effect of impurities for the 5 = 1/2 Heisenberg antiferromagnetic chain has been studied 
extensively. Eggert and Affleck [|l],|[ have studied the open chain, where edges cause the 
inhomogeneity. They investigated the static structure function and the local susceptibility, 
using the conformal field theory and a Monte Carlo method. Laukamp, et al. have also stud- 
ied the open chain problem by a Monte Carlo method and a DMRG method 0. Magnetic 
structures at edges cause interesting effects on the NMR line shape 0,|]. In particular, the 
observation 0] of a broad background in the NMR spectrum in Sr2Cu03 has been found to 
agree with the features predicted by the theoretical work on the field-induced local staggered 
magnetization p[. The temperature dependence of the susceptibility of an ensemble of open 
finite chains with various lengths has been also studied experimentally and theoretically [Q . 

The properties for 5 = 1/2 chains with random exchange coupling have been also studied. 
Random distribution of the interaction brings a new type of low temperature phase such 
as random singlet phases |^ and furthermore various other types of randomness- induced 
phases have come out according to the distribution of the bonds [§,0. 

Impurity effect in the S = 1 AF chain has been also investigated. The existence of the 
edge state is one of the most interesting properties of the Haldane state [0 which was 



exactly pointed out in the AKLT |TI| chain. This edge state has been studied in detail 
by numerical method |1^,|13|, where the singlet and triplet (Kennedy triplet) states become 
degenerate in the thermodynamic limit (the four-fold degeneracy). Furthermore a doping of 
an impurity of S* = 1/2 brings a local magnetic structure which causes an interesting energy 
structures [|14| -[16[|. The effects of bond impurity has been also studied [|17|. There is the 



robustness of the local structure due to the 5* = 1/2 impurity, while the interaction between 
the local structures has been pointed out to be very weak because of the quasi degenerate 
energy structure [|1^ 



In this paper we study magnetic structures in S" = 1/2 chains due to various spatial 
configurations of one or two bond impurities in the uniform system. Here we consider a 
uniform antiferromagnetic Heisenberg chain and we regard a weaker or stronger bond than 
the bulk bonds as an impurity bond. First we review the properties of the uniform antifer- 
romagnetic Heisenberg model in open chains and investigate the temperature dependence 
in detail. When we put an impurity in an open uniform chain, it is found that such an 
impurity effectively separates the system and the right domain and the left domain behave 
almost independently at modestly low temperatures. According to whether the number of 
spins of the domain is even or odd, magnetic property of the domain shows very different 
characteristics. We investigate the magnetic properties by calculating the magnetization 
profile, the spin correlation function, the local susceptibility and the local field susceptibility 
(see next section). 

We also study bond-impurity effects in the bond-alternating chain, where 'impurity' 
means a defect of the alternating order. In the bond-alternating system, the correlation 
length is finite and similar to the case oi S = 1 Haldane systems. The impurity induces a 
localized magnetic structure around it. 

In some cases the position of the bond impurity can move. We study the force acting 



between the impurity bonds. We find the force is attractive. 

This paper is organized as follows. In the next section, we explain briefly the method 
used in this study. In Sect. |T|, we investigate magnetic structures for the uniform systems 
with impurities. In Sect. |^, the effects of bond impurities for the magnetic structures of 
bond-alternating chains are studied. In Sect. |V], we study the force between bond impurities. 



Sect. VI is devoted to the summary and discussion. 



II. MODEL AND METHOD 

In this paper we study the low temperature properties of 5 = 1/2 antiferromagnetic 
Heisenberg chains with bond impurities. The Hamiltonian is generally given as 

i i 

where Sj = (S'f , Sf, S^) are the S = 1/2 spin operators. In Sect III, we study the uniform 
chain where Ji = J except at the position of the impurity bonds, j, where Jj is different 
from the bulk value J. Hereafter we take J as the unit of the energy. In Sect IV, we study 
the bond-alternating chain where Jj takes a strong bond Ji and a weak bond J2 alternately. 
There the impurity means defects of the order of the alternation. 

We mainly use the loop algorithm with continuous time quantum Monte Carlo method 
(LCQMC) [J19l in this study. This method overcomes the problem of long autocorrelation 
in Monte Carlo update. Furthermore replacing discrete time with finite Trotter number, an 
algorithm using continuous time has been introduced [^ and the nuisance of the extrap- 



olation of the Trotter number has been released. These improvements allow us to study 
systems at very low temperatures [^. In the present work, we performed 10^ Monte Carlo 
steps (MCS). Here a MCS means a update of whole spins. 

In the World Line quantum Monte Carlo method (WLQMC) we can control the value of 
total magnetization [Mz) in the initial state and keep it by suppressing the global flip which 
changes the magnetization. Thus, generally it is not difficult to obtain the true ground state 
configurations fixing the value of the total magnetization. However, in a standard LCQMC 
the number of world line and the winding number are updated automatically and we can not 
specify the value of the magnetization M^. We adopt LCQMC with fixed total Mz values. 
We perform the standard LCQMC and store the data separately according to M^. When 
we need information for a specific value of Mz, we use the data with that value of Mz only. 
This method worked successfully in the study of the site impurity problem {S = 1/2 spin) 



in the 5* = 1 Heisenberg antiferromagnetic chain [|18|. In order to check the method in the 
case oi S = 1/2 chains, we confirmed the agreement between the low temperature results 
obtained by the LCQMC for small sizes such as L = 20 and 21 and those obtained by exact 
diagonalization. 

We study the local magnetic properties by investigating the following properties: the 
local magnetization 

m. = (S!), (2) 

and the spin correlation function 



C{t,j) = {StS]). (3) 

At finite temperatures, no local magnetization appears, rrii = 0. For the purpose of 
detecting the local magnetic structure, we calculate the local susceptibility introduced by 
Eggert and Affleck: 



X. - |(5f) 



= l3j2{S^Sr}, (4) 

h=Q j 



where h is the uniform field {hi = h) and (3 = 1/T. This quantity is nonzero even at 
finite temperatures where all M^ subspaces contribute. For the odd-chain (a chain with odd 
number of spins) Xi is proportional to the local magnetization mj at T = 0. 

We also study another type of local susceptibility, which we call local field susceptibility 
in this paper. 
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\st{r)St)dT - P{S-t)\ (5) 

hi=Q 



which has been investigated as an indication of quantum fiuctuation in the study of quantum 



spin glasses or diluted models |22 



III. BOND IMPURITIES IN UNIFORM OPEN CHAINS 

In this section we study the effect of bond impurities in the uniform S = 1/2 chain. 
First, we review the magnetic properties of the open chains and investigate the temperature 
dependence of them. Next, we investigate magnetic properties in open odd-chains with one 
bond impurity. And then, we investigate systems with two impurities and how magnetic 
properties are changed by shifting the position of the two impurities. 

A. Open uniform chain 

Eggert and Affleck [||,|| studied characteristics of magnetic properties for open uniform 
even-chains at low temperatures. They found that local susceptibility remains near the 
edges for even-chains at modestly low temperatures. Laukamp, et al. recently studied 
the magnetic properties for odd- or even-chains in a fixed M^ space by DMRG. They found 
that the local susceptibility forms a simple sinusoidal form. These features are relevant to 
study the temperature dependence of experimental results. 

Fig. ^(a) and (b) show Xi ^^ ^ very low temperature T = 0.01 of chains with L = 63 
and 62, respectively. The odd-chain shows a nodeless shape of Xi while Xi of the even-chain 
shows a node, which reproduce the result of Laukamp, et al. As the temperature increases, 
Xi at middle of the chain shrinks as shown in Fig. ^ (L = 63, T = 0.05). In a longer chain, 
the shape reproduces that of Eggert and Affleck [|I|,g ( L = 128, T = 0.067, Fig. |(a)). 
At this temperature Xi of ^in odd-chain [L = 127) shows a very similar shape to that of 
L = 128 as shown in Fig. ^ (b). On the other hand at the very low temperature T = 0.01, 
the profiles of {xi} for odd-chain and even-chain are very different as shown in Fig. |] (c) 
and (d). These profiles are essentially the same as those in Figs. ^ 



Here it should be noted that the amphtude of Fig. |T](b) is very small. In Fig. ^ we 
show the data of x« at T = 0.05 where we find almost the same shape but the amplitude is 
about 20 times larger than that of T = 0.01. This temperature dependence is understood as 
follows. In the subspace of the total magnetization M^(= J2i ^^i) = 0, the local susceptibility 
vanishes by the definition Eq. (^. In the ground state we expect that Xi=0 assuming the 
ground state is singlet. If a triplet state exists above it (with the energy gap AE) and 
the other states have high energies and do not contribute to the thermal distribution, the 
amplitude of the shape is given by 

X. = 7X7(M, = 1), (6) 

where 

2g-/3A£; 

^ = l + 3e-/3AS' ^^'> 

and xii^z = 1) is the local susceptibility in the lowest level of M^ = 1. We suppose that Xi 
is proportional to that in Fig. 0(b). 

This ratio 7 is obtained by the fraction of MCS in the subspaces of M^ = ±1. At 
T = 0.01 the fraction is about 0.34 % and all other MCS are in the subspace of M^ = 0. 
Thus the above assumption that only the lowest singlet state and triplet states contribute 
to the thermal distribution is confirmed to be valid. 

At T = 0.05 the fraction of Mz = ±1 is observed to be about 34 % in MC simulation 
and small number of MCS distribute in the spaces of \Mz\ > 2. Thus the distribution is 
about 100 times larger than that at T = 0.01, which is consistent with the ratio, Xi{T = 
0.05)/x*(T = 0.01) ~ 20, taking into account /3 in the definition of Xi (Eq. (|)). 

Using the relation Eq. (^, from this fraction 7 = 0.0034 at T = 0.01 the energy gap is 
estimated as 

AE{L = 62) ~ 0.064. (8) 

In the chain of L = 128 the energy gap is similarly estimated as 

AE{L = 128) ~ 0.031. (9) 

The energy gap becomes small as the size increases proportionally to 1/L as is expected for 
S = 1/2 antiferromagnetic Heisenberg chains. 

In Fig. |^(a) and (b), x'f'^"''' fo^' -^ = 63 and 62 are shown, respectively. In both cases, 
x'°™' also shows a zigzag shape, although in Fig. |^(b) there is a node as in Fig. |l] (b). This 
zigzag behavior indicates that the amplitude of quantum fluctuations is larger at sites where 
the local spin is antiparallel to the field. 

B. One bond impurity 

As bond impurities we consider three types of inhomogeneity: (a) a strong bond impurity 
(Js=2), (b) a weak bond impurity (J„=0.5), and (c) a complex structure of impurities 
consisting of both strong (Js=2) and week (7^=0.5) bond impurities. Exactly speaking, the 



last type does not consist of a single bond impurity but we investigate this case here because 
of an analog to the formers. In Fig. ^(a), (b) and (c), the magnetic profiles {rrii} obtained 
by LCQMC at T = 0.01 in the M^ = 1/2 space are shown for the types. Configurations 
of bonds are also illustrated in the figures. These profiles are considered to represent well 
the ground state because the structure shows little dependence on the temperature within 
a fixed value of magnetization (M^). In Fig. |^(a) the left hand side of the impurity shows a 
steady staggered structure, while there remains little magnetization in the right hand side. 



This structure is easily understood from the viewpoint of the Lieb-Mattis theorem ^3 
S = 1/2 open odd-chain has the doublet M^ = ±1/2 state as the ground state and the local 
magnetization is nonzero, while the even-chain has the singlet M^ = ground state and no 
local magnetization appears. In the configuration of Fig. ^(a), the spins connected by the 
strong bond impurity form a singlet state and the whole system is efficiently divided into 
two parts consisting of spins lst-29th (odd) and 32nd-61st (even). Thus, approximately, the 
doublet ground state appears in the left domain and the singlet ground state appears in the 
right domain. On the other hand in the configuration of Fig. |^(b), the weak bond itself 
cuts the whole chain and the system is divided into two part: lst-30th (even) and 31st-61st 
(odd). Consequently the structures of the right hand side and the left hand side exchange. 
Due to the same reason, a magnetic structure also appears in (c), where the two domains 
are more clearly separated and the contrast between the right and left sides is more clear. 

In order to investigate the correlations in each domain, we calculate the two-point cor- 
relation functions of the model of Fig. ^ (a) in the fixed M^ = 1/2 space at T = 0.01. 
The two-point correlation functions from the left edge C(l, i) and the right edge C{i, L) are 
compared in Fig. ^(a) and the two-point correlation functions from the middle site of the 
left domain C{15,i) and the right domain C(i,47) are compared in Fig. |^(b). While the 
magnetic profile {rrii} of each domain in the ground state is very different from each other, 
almost the same correlation in both domains exists which is also very similar to that in the 
pure S = 1/2 chain. Because the correlations are almost the same, we expect that the spins 
in the right domain fluctuate coherently keeping the mutual correlation. That is to say, the 
spins in the right domain have a similar magnetic proflle to that in the left domain in some 
direction and the direction is fluctuating. The same tendency is also observed in the models 
of Fig. I (b) and (c). 

We also investigate Xi and Xi"™' of the lattice with one strong bond impurity for L = 127, 
where the impurity bond locates between 64th and 65th sites. Here the left and right 
domains contain 63 and 62 sites, respectively. When the temperature is modestly low such 
as T = 0.05 we flnd that each domain independently shows magnetic property according to 
whether the number of spins is odd or even, as shown in Fig. |^. At a very low temperature 
such as T = 0.01, {xi} in the right hand side is expected to be given by that in Fig [l| 
(b) where the amplitude is very small with a node, if the domain behaves independently. 
On the contrary, the amplitude of Xi in the right hand side is enhanced and the one-node 
structure disappears, as shown in Fig. ^. This change is understood as follows. If the domain 
would behave completely independently, Xi is nearly zero because the singlet ground state 
dominates and the amplitude of the one-node structure decreases as Eq. @ as we see in the 
previous subsection. However in the present lattice, the right domain can be regarded as an 
even-chain only approximately because of nonzero interaction through the impurity bond. 
Thus the magnetization of the domain is no more a good quantum number. Consequently in 



the lowest state of the right domain, Xi may have some structure. The structure is expected 
to have small local magnetic moment because the state would be nearly singlet. In Fig. |l^, 
we show Xi cit an intermediate temperature. Although we do not show it here, the magnetic 
profile {rrii} is found to be proportional to {xi} in this lattice. 

The local field susceptibility at T = 0.05 is shown in Fig. ^ (b). When we decrease 
the temperature to T = 0.01, {xi"'^"'''} is found to change smoothly to a shape given by the 
combination of Fig. ^ (a) and (b) in each domain. 

C. two bond impurities 

In this subsection we consider two-impurity problem. There are various configurations 
with two bond impurities concerning to combinations of types of impurities and also to 
whether the number of spins in each domain is odd or even. Here we deal with two types of 
lattices of L = 61 [J = 1) with two strong bond impurities (Jg = 2). The impurities divide 
the lattice into three domains: model (A) a lattice where the domains contain successively 
odd, even, and even number of spins and model (B) a lattice where the domains contain odd, 
odd, and odd number of spins. The lattice structures and the magnetic profiles at T = 0.01 
are presented in Fig. |ll](a) and (b). These magnetization profiles represent well the ground 
states of the models because almost all (100 and 98.9%) of Monte Carlo steps are distributed 
in the Mz = ±1/2 space which means the system is in a doublet state. The feature of these 
magnetic profiles can be understood from the analogy to the arguments given in the previous 
section. Approximately, the doublet ground state appears in domains with odd number of 
spins, while the singlet ground state appears in domains with even number of spins. 

For the model (B), in a naive picture we may consider the magnetizations of domains to 
be +1/2, —1/2 and +1/2. Hereafter we denote this configuration by (1/2,-1/2,1/2). Fig. 



12| shows the summation of magnetization per site from the left edge site 

Mz{j)=J2mi, (10) 

which shows that the right and left domains have positive magnetizations, while the middle 
domain has a negative magnetization. In this figure we find tendency of a staggered domain 
magnetization. The magnetizations of the domains are reduced to 0.3, —0.1, and 0.3, respec- 
tively. This reduction comes from a linear combination of the state (1/2, 1/2, —1/2) and 
(—1/2, 1/2, 1/2) due to the quantum fluctuation. Looking on the domain magnetization 
as an effective 5* = 1/2 spin interacting by an exchange J, this system is modeled by a 
three-site Heisenberg model H = JSi ■ S2 + JS2 -83. In the Mz = 1/2 space the eigenvalues 
and the eigenvectors are described as 

El = -J, 

\h) = i/VqH + +-) - 2| + -+) + I + +^)) 

^2 = 0, 

102) = i/v^(| + +-)-!-++)) (11) 

E3 = J/2, 

103) = 1/V3(| + +-) + I + -+) + I - ++)). 



Here, we find (0i|S'f |0i) = 1/3 and (0i|S'||0i) = —1/6. The reduced magnetization of the 
domains is close to this set (1/3, —1/6, 1/3) and thus the domain structure is well described 
by the state 0i. We find that the state of each domain can be regarded as a doublet state 
approximately. 

Finally in Fig. |T3| we show Xi at modestly low temperature (T = 0.05) for a long chain 
of the model (A), where we find again the independent behavior of the domains. 

IV. BOND IMPURITY IN BOND- ALTERNATING CHAINS 

In this section, we investigate models with bond alternation, ■ ■ ■ JiJ2JiJ2- ■ ■, where 
Ji > J2. Here we study the effect of a defect of the alternation, such as 
■ ■ ■ J I J 2 J I J 2 Ji J I J 2 Ji J 2 J I ■ ■ ■• As mentioned in the introduction, the alternating chain has 
a energy gap between the ground state and the first excited state and the spin correlation 
length is finite. These features are similar to those in the S = 1 AF Heisenberg model which 
is the Haldane system. As well as in the case of the Haldane system, there are edge states in 
this model, too. The edge state strongly depends on the bond situation at the edge. Thus 
we study the following four cases of the configurations of the bonds at the edges and the 
center: 

(a) the two strong bonds are at the center and the edges terminate with strong bonds, 

(b) the two strong bonds are at the center and the edges terminate with weak bonds, 

(c) the two weak bonds are at the center and the edges terminate with weak bonds, 
and 

(d) the two weak bonds are at the center and the edges terminate with strong bonds. 
Here we take the strong bond to be Ji = 1.3 and the weak bond to be J2 = 0.7. The mag- 
netic profiles of (a)-(d) are drawn in Figs. ^^ where simulations were performed at T = 0.01 



and in the M^ = 1/2 space. A magnetization is induced locally around the impurity. 

First we consider the cases where the bonds at edges are strong, i.e., model (a) and (d). 
If we allocate a singlet pair at each strong bond, a structure with neighboring two strong 
bonds remains at the center in the model (a), while one site remains in the model (d). The 
magnetization of M^ = 1/2 is assigned in the remaining part of the lattices to induce a local 
magnetic structure. Because the edge bonds are strong, no magnetization is induced at the 
edges. The ground state is simple doublet and 100% of Monte Carlo steps are distributed in 



the Mz = ±1/2 space. Fig. 14(a) and (d) are considered to describe well the magnetization 
profiles of the ground state of both models. In the case (a) negative magnetization appears 
at the middle site, while in the case (d) positive magnetization appears there. This structure 
is naturally understood as follows. The interaction of the three spins at the center of the 
model (a) is represented as three-spin model coupled by the strong bonds where the state 
|0i) in Eq. ( ITll) gives the ground state, while in the model (d) a spin at the center is isolated 
from the others which form singlet states. 

On the other hand, in the model (b) and (c), if we allocate singlet spins at the strong 
bond, spins remain at the center and also at both edges and thus there are three positions 
for magnetic moments. Indeed in both (b) and (c) models, magnetization is induced locally 
around the impurity and the right and the left edge sites. Fig. |T3| shows the summation 



of magnetization per site from the left edge site for the model (b). In Fig. |T5| the values 



of the left plateau and right plateau are 0.166 and 0.333, respectively. From this figure we 



find a spin 1/6 locates at each local structure. This deceptive fractional magnetization is 
considered to come from mixing of states. As we mentioned in Sect. Ill, the eigenstates in 
Mz = 1/2 of the open chain of three spins are given by Eq. (|11]), where 

1/3 and 

and (12) 
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(031^3103) = 1/6. 

Although there are energy gaps between the state |0i), |02), and |03), the tempera- 
ture is considered to be much larger than these energy gaps and these states are almost 
degenerate and thus the expectation values of S'f and 5*1 are (l/3+0+l/6)/3=l/6 and 
(-l/6+l/2+l/6)/3=l/6, respectively. 

In order to check that the ground state of the type (b) is represented by the |0i) of the 
three-spin model, we investigate a short chain of L = 21 by the diagonalization method. 
Because the length of the chain is short, we choose a shorter localization length. For this 
purpose, we set strong and weak bonds Ji = 2 and J2 = 0.5, respectively. Fig. ^ (a) 
shows the magnetic structure in the ground state of this model and Fig. |T^ (b) shows the 
summation of the magnetization per site from the left edge site. The net magnetization 
around the right and the left edge is positive and that around the impurity site is negative. 
The value for the left plateau of the Fig. [l^(b) is about 1/3 and this corresponds to 
((/'ilS'fl^i), and the value of the second plateau is about 1/6, which also corresponds to 
{(l)i\Sl\(t)i) + (^il^ll^i). Thus the three-spin model is valid for this model of L = 21. On 
the other hand in the case of L = 63, the energy gap is much smaller than the temperature. 
Thus all the quasi-degenerate states contribute to the observation. 

This degeneracy is also seen in the distribution of M^. M^ distributes as shown in 
Fig. O. Here the distribution for M^ = 3/2 is about 12.5 %. There are three states in 



the Mz = 1/2 space and one state in the M^ = 3/2 space in the three-spin model. The 
fact that the distribution for total M^ = 3/2 is about 12.5 % means that T = 0.01 is much 
higher compared with the energy gaps between these four states. In the three-spin model the 
eigenvalue of the state of M^ = 3/2 is J/2, while the eigenvalues of the state of M^ = 1/2 are 
—J, 0, and J/2 (Eq. (pA])). In principle we can obtain the energy gap from the temperature 
dependence of the distribution. However it is too small to be detected here. Thus these three 
states in the M^ = 1/2 space appear in equal probability, and (S^) = (^f) = {S^) = 1/6 in 
the Mz = 1/2 space as mentioned above. In the M^ = 3/2 space we observe S* = 1/2 moment 
at each local structure as shown in Fig. |T^(a) and (b). We find a similar scenario for the 
model (c). Thus we conclude that in rather strong bond-alternating systems, local magnetic 
structures induced by a bond impurity or weak edge bonds have an effective 5* = 1/2 spin 
and they behave almost independently. 

V. FORCE BETWEEN BOND IMPURITIES 

If we allow the position of impurities to move, impurities move by the force between 
them. In the equilibrium state, the distance between impurities distributes in the canonical 

9 



distribution for the interaction energy between the impurities. Here we study the force 
between the bond impurities. 

A. Two bond impurities in the uniform chain 

In this section we consider the force between bond impurities in the uniform chain. First, 
we deal with L = 12 periodic chains (J = 1) including two strong bond impurities (Js = 2) 
at various positions. We study the ground state energy of the system as a function of the 
position of impurities. An impurity is put on the first bond and the other is put on the n-th 
bond. We calculate energies of the system as a function of the distance 

A = n-1, (13) 

which is shown in Fig. |T^. The energy -E'(A) corresponds to the potential energy of the 
force between the impurity bonds. We find that the cases in which A is even are more 
energetically favorable than those of odd A. Among the cases of even-distance the energy is 
lowest when the distance is shortest, which indicates that an attractive force acts between 
the bond impurities. 

To confirm this conclusion for the attractive force in a larger system, we preformed 
Monte Carlo simulations for an L = 60 periodic chain. Here we use the following algorithm 
to update the system. We start with an arbitrary configuration with two impurities, and 
then we update the spin variable with the standard LCQMC Next we update the bond 
configuration in the following way. We exchange the neighboring bonds sequentially from 
the left. If the two bonds are the same, we skip to the next. We choose a new configuration 
of bonds by the thermal bath algorithm, i.e., we take the exchanged configuration with the 
probability: 

W 
P = — (14) 

where W is the Boltzmann weight for the original bond configuration and W is that 
of the exchanged configuration for the given spin configuration \a): W={a\e^'^''^^\a) and 
iy=(cr|e~^^B'|o"). We plot the distribution of the distance between the impurity bonds in 
Fig. ^. We expect that the distribution should be the canonical distribution with the po- 
tential energy of the interaction of the impurities. We see in Fig. ^ that the most probable 
distance is A = 2. Thus we conclude that A = 2 gives the lowest energy. When A is 
odd, the probability is low, which indicates that the energies for odd-distances are higher 
than those for even- distances. Among the even As, the probability decreases as the distance 
becomes long, which indicates the attractive force acts. 

We also studied the case of weak impurity bond and found that the attractive force acts 
similarly to the case of the strong impurity bonds. 



B. Two defects in an alternating chain 

In the alternate chain (■ ■ • Ji J2^i^2 ■ ■ ■ )' ^^e system may have a pair of defects by shifting 
a position of a strong bond by one. 

10 



■■■JlJ2 JlJlJ2J2 JlJ2---. (15) 

If we shift the position furthermore, the system has a configuration 

■ ■ ■ J1J1J2J1J2J2" ' ■> (16) 

etc. We study dependence of the energy on the distance (A) between the positions of Ji Ji 
and J2-^2- In the case of no defect we define A = 0, A = 1 for Eq. (p!5|), A = 2 for Eq. (|16D 
and so on. In Fig. |2T] we plot the ground state energy as a function of A obtained by exact 
diagonalization for L = 24 with Ji = 2 and J2 = 1- Here we find again an attractive force 
between the impurities. 

VI. SUMMARY AND DISCUSSION 

In this paper we explored effects of various bond impurities on the low temperature mag- 
netic properties in the uniform and bond-alternating antiferromagnetic Heisenberg chains. 

In particular we clarified the temperature dependence of the magnetic structure of the 
open uniform chain. At modestly low temperatures systems show a similar local suscepti- 
bility profile regardless of whether the number of spins is even or odd. On the other hand 
at very low temperatures, the profile is very different from each other. 

We also studied domain structures which are separated by bond impurities in the uniform 
chain. There we found that each domain behaves almost independently at modestly low 
temperatures. However, at very low temperatures in the domain with even number of spins, 
a magnetic structure appears due to the interaction between the domains, although such 
structure does not exist in a completely isolated open chain with even number of spins. Thus 
the interaction between the domains is not negligible at very low temperatures. We also 
found that the interactions between domains with odd number of spins can be approximately 
modeled by a simple system with S" = 1/2 spins representing the domain magnetization. 
The effective interaction J of the model is not necessarily small at T = 0.01. Because the 
energy gap due to the finite size in the uniform chain is proportional to l/L, we expect that 
J also becomes small as 1/L. 

Effect of bond impurities are investigated also in the bond-alternating chain. The role 
of the bond impurity is very different from that in the uniform chain, i.e. bond impurities 
induce a magnetic structure around them, while they cause the division into domains in 
the uniform chain. In the bond-alternating chain the effect is quite similar to that of the 
site impurities of 5* = 1/2 in the S = 1 antiferromagnetic chain because of the gapful 
nature. Here induction of magnetic moments at edges depends on whether the edge bonds 
are strong or weak. In the bond-alternating chain, the magnetic structures behave almost 
independently even at very low temperatures, contrary to the case of the uniform chain. If 
we model the interaction of the magnetic structure by a effective spin model, the interaction 
J is expected to be exponentially small with the system size, which is due to the gapful 
nature of the bond-alternating chain. 

We also investigated the force acting between the bond impurities and studied the dis- 
tribution of impurities when impurities are allowed to move. It turned out that the force is 
attractive both between the impurity bonds in the uniform chain, and between the defects 
of alternation in the bond-alternating chain. 

11 



We hope that the present study helps to analyze the magnetic properties at low temper- 
atures such as observed by NMR measurement. 
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FIGURES 
FIG. 1. Local susceptibility Xi at T = 0.01. (a) L = 63 and (b) L = 62. 

FIG. 2. Local susceptibility Xi at T = 0.05 for L = 63. 

FIG. 3. Local susceptibility Xi at T = 0.067 (a) L = 128 (b) L = 127 and at T = 0.01 (c) 
L = 127 and (d) L = 128. 

FIG. 4. Local susceptibility Xi at T = 0.05 for L = 62. 

FIG. 5. Local field susceptibility x^°'"'^ at T = 0.01 for (a) L = 63 and (b) L = 62. 

FIG. 6. Local magnetization in the M^ = 1/2 space at T = 0.01 for the models with L = 63 
(a)-(c): (a) strong bond impurity (Jg = 2) (b) weak bond impurity (Jw = 0.5) (c) strong (Jg = 2) 
and weak (J^ = 0.5) bond impurities. 

FIG. 7. (a) Two-point correlation functions from the left edge and the right edge are compared, 
(b) Two-point correlation functions from the middle site of the left domain and the right domain 
are compared. 

FIG. 8. (a) Local susceptibility of the lattice with one strong bond impurity for L = 127 at 
T = 0.05, where the impurity bond locates between 64th and 65th site, (b) Local field susceptibility 
of the lattice with one strong bond impurity for L = 127 at T = 0.05. 

FIG. 9. Local susceptibility of the lattice with one strong bond impurity for L = 127 at 

r = o.oi. 

FIG. 10. Local susceptibility of the lattice with one strong bond impurity for L = 127 at 
T = 0.02. 

FIG. 11. (a) Lattice structure and local magnetization for the model (A), (b) Lattice structure 
and Local magnetization for the model (B). 

FIG. 12. Summation of magnetization per site from the left edge site for the model (B). 

FIG. 13. Xi at T = 0.05 for a long chain (L = 127) of the type of the model (A). 
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FIG. 14. Local magnetization of the models (a)-(d) at T = 0.01 in the Mz = 1/2 space, (a) 
and (c) contain 63 sites and (b) and (d) contain 65 sites. The diamonds denote the strength of 
bonds Jj, those at high positions denote Ji and those at the low positions denote J2. Details are 
shown in text. 



FIG. 15. Summation of magnetization per site from the left edge site for the model (b). 

FIG. 16. (a) Local magnetization of a bond-alternating chain with L = 21, which is the same 
type as that in Fig. |l^. (b). The strong and weak bonds are Ji = 2 and J2 = 0.5, respectively, 
(b) Summation of the total magnetization per site from the left edge. 

FIG. 17. Distribution of M^ for the model (b) in Fig. ||. 



FIG. 18. (a) Local magnetization for the model (b) in Fig. 14 in the M^ = 3/2 space at 



T = 0.01. (b) Summation of the total magnetization per site from the left edge. 

FIG. 19. Ground state energy as a function of A obtained by exact diagonalization for L = 12 
periodic chains with two strong bond impurities. 

FIG. 20. Distribution of A by a Monte Carlo method for L = 60. Open squares and circles 
show data for T = 0.1 and T = 0.05, respectively. 

FIG. 21. Ground state energy as a function of A obtained by exact diagonalization in L = 24 
for bond-alternating systems with defects. 
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